Preserving universal resources for one-way quantum computing 
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The common spin Hamiltonians such as the Ising, XY, or Heisenberg model do not have ground 
states that are the graph states needed in measurement-based quantum computation. Various 
highly-entangled many-body states have been suggested as a universal resource for this type of 
computation, however, it is not easy to preserve these states in solid-state systems due to their 
short coherence times. Here we propose a scheme for generating a Hamiltonian that has a cluster 
state as ground state. Our approach employs a series of pulse sequences inspired by established NMR 
techniques and holds promise for applications in many areas of quantum information processing. 
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Measurement-based quantum computation (MQC) is a 
new computing paradigm [lj. Of particular interest are 
universal resources of one-way quantum computation, a 
MQC scheme that requires only local measurements [2j. 
In the original scheme of one-way quantum computing, 
one initially creates a many-qubit cluster state by apply- 
ing phase gates or equivalent gate operations which can 
be realized using the Ising interaction between qubits. 
Many promising methods to generate cluster states us- 
ing solid-state qubits have been proposed [3, 4]. However, 
these states are not the ground states of spin Hamilto- 
nians with typical qubit-qubit interactions of Ising, XY, 
and Heisenberg types [5|- 

One of the established universal resources are two- 
dimensional (2D) cluster states. Another promising can- 
didates is the Affleck-Kennedy-Lieb-Tasaki (AKLT) state 
on the honeycomb lattice a resonance valence bond 
(RVB) type state which is a special projected entan- 
gled pair state (PEPS) @, i]. Yet, the AKLT state re- 
quires non-trivial Hamiltonians with spin greater than 
1/2, which are not easy to realize in solid-state systems. 

In this Letter, we present a new method of constructing 
stabilizer Hamiltonians whose ground states are the de- 
sired cluster states. Our approach relies on manipulating 
a two-body Hamiltonian using pulse-sequence techniques 
developed in the nuclear magnetic resonance (NMR) con- 
text [§, We show that, starting from the Ising and 
XY models, one can induce an effective dynamics de- 
scribed by a stabilizer Hamiltonian [2| 
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(i) 



where K% = erf ®j en bhd(i) a j are the correlation opera- 
tors and the direct product runs over all nearest neigh- 
bors of the lattice site i (erf and erf are the Pauli ma- 
trices). Thus, the main idea of this paper is to gener- 
ate a stabilizer Hamiltonian by transforming the original 
Hamiltonian. 

We assume the original Hamiltonian to be of the form 



H = Hq + Hi n t where 



(2) 



is a single-qubit part and H m t the interaction part. We 
take H int to be of Ising ffi sing = Y*<j J ij a i a j > XY 
Hxy = J2i<j Jij[XY]ij, and Heisenberg form H H = 
J2i<j Jij[XYZ]ij. Here, we set K = 1, and use the 
shorthands [JCY]y 



<T?0-j + <jf(Tj 



a fa* 



a\a) and [XYZ] t 



afa*. We will assume = J if i and j 



are nearest neighbors and Jij = otherwise. 

Note that a single correlation operator can be obtained 
using a single-qubit Hamiltonian. For example in a one- 
dimensional (ID) qubit array, K2 = erf erf erf can be gen- 
erated by the time evolution operator e^/ 2 * 1 ^ 1 +a2 +a3 ' . 
However, it is not evident how to obtain a sum like 
K2 + K3 from the single-qubit Hamiltonian. 

Most fabricated solid-state qubit systems are nano- 
devices, because a smaller size makes them more ro- 
bust to decoherence. An example are quantum dot sys- 
tems where smaller dots have larger energy-level spac- 
ing. Since with diminishing size it becomes difficult to 
address these devices individually, it is of interest to con- 
sider switching on/off Hq and H mt independently. We 
will show that, by using appropriate pulse sequences, this 
is possible even if we start from an always-on Hamilto- 
nian [m . 

Ising model- We now show how to construct the sta- 
bilizer Hamiltonian using the relation 



'- ale"' 

2/ iflcr'crj 



cos(26»K +sin(20)aj'CTf, 
cos(20)erf - sin(20)oi erf . (3) 



An important consequence of these equations is that, 
for 9 = 7r/4, we can increase the order of the Pauli- 
matrix terms as in e~ l T <T i cr 2 C r^ e l f a i a 2 



erf "2 and 



-erf erf- F° r a ID -/V-qubit chain, 
the starting single-qubit Hamiltonian is given by 

H s = n(af+ >>' - rr'(i. (4i 
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By applying H int as 

— iOy^ -, (7 z (j z tt i& Y" 1 ■ - ■ erf <rf 



(5) 



we obtain the ID stabilizer Hamiltonian -ff s tab for & — 
7r/4. As an example, for N — 3 qubits, starting from 
i/ s = Q,(a\ + erf + cr|) we obtain i7 sta b = Q(-o~fcr% — 
cr i f7 2 f7 3 ~ <7 3 a 2)- 11 the system of A qubits has peri- 
odic boundary conditions, we start from the Hamilto- 
nian H s = ft Y]^-i &f- Since Ising-type interaction terms 
commute and the time-evolution operator in ([5]) factor- 
izes, this process can be straightforwardly extended to 
2D and 3D qubit systems, thus realizing the universal 
resource discussed in the introduction. Consequently, for 
Ising interactions, we can construct the stabilizer Hamil- 
tonian by switching on -ff ising only once. 

We will now implement the above equations into a 
pulse sequence that is familiar in the NMR context 
We assume that each pulse is sufficiently strong such 
that interactions between qubits can be neglected dur- 
ing the pulse sequences. We describe the time evolution 
of the system by the density operator pit) whose time de- 
pendence is given by p(t) = exp(—iHt)p(0) exp(iiit) for 
time-independent H . It is convenient to use the follow- 
ing schematic notation for this evolution: p(0) p(t). 
Then the process 



P(P) 



TlH U 



-TtHu 



Pit) 



rH B , 



(6) 



for T\ = 7r/(4J) corresponds to p(0) p(t) where 

t = t + tt/ (2 J) and r can be chosen arbitrarily. Note that 
at the physical time t the state of the system is obtained 
from the initial one by the time-evolution operator 



-irH s , 



/4) E, 



' e 



(7) 

As illustrated by Fig. Ufa), H sta b becomes the effective 
system Hamiltonian. Its ground state is the originally 
prepared cluster state, which is therefore preserved. 

XY model- Next, we show how to generate the sta- 
bilizer Hamiltonian using the XY interaction, assuming 
that H and H mt can be switched on/off independently. 

The stabilizer Hamiltonian is formed step by step 
by bonding the nearest-neighbor operators. This 
is because the XY interactions do not commute, 
\[XY\i-i t i, [Xy] iii+1 ] ^ 0. We start from 



e -ie[XY] 12a x e ie[XY] 12 = 

e -i0[XY] 12(T y e WlXY] 12 = 
e -i9[XY] 12(T z e ie[XY]i 2 = 



cos(20)of - sin(20)ofaf, 
cos(20)of + sin(20)o-f of , 
cos 2 (20)crf + sin 2 (20)erf 

i S m(40)[aM-^f]- (8) 



For 9 = 7r/4, these transformations increase the order of 
the Pauli-matrix terms as o~f — > and a\ — ¥ o~\o~ 2 . 

For a\ one obtains <r\ — > erf. 
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FIG. 1: (a) Switching on/off the parts H 3 and Hi n t of the 
Hamiltonian H = H s + Hint gives rise to an effective time 
evolution described by the stabilizer Hamiltonian i/ sta b [see 
Eq. (0] whose ground state is a cluster state, (b) Graph 
representation of a one-dimensional 4-qubit (squares) cluster 
state stabilized by the Ising interaction, (c) Two-dimensional 
4x4 twisted cluster state stabilized by the XY interaction, 
a universal resource for one-way quantum computation. The 
dashed lines and Pauli operators in each direction illustrate 
the twistedness of the state and the corresponding stabilizer 
Hamiltonian. 



We now show how to construct a 2D stabilizer Hamil- 
tonian. First we construct the ID stabilizer Hamiltonian, 
starting from 



H s = Q(-af + a\ 



E 

i=3.JV-2 



(9) 



In the specific case of six qubits in ID, by applying Eq. 
to [XY} 12 , [XY] M , and [XY] 56 , we obtain: 



-iSi 



H s e 



= n(ffj a\ + al erf + crfcrf + crfcrf 

+ alai + alat), 



(10) 



where S\ = j^2i=i 3 [XY]{ 2 i-i,2i}- Repeating this step 
with S 2 = f J2i=iA XY ]{2i,2i+i}, we get the ID sta- 
bilizer Hamiltonian H m = e~ iS2 e~ iSl H s e iSl e iS2 that 
reads explicitly 

U Ol^r Z ^ Z ^r X I ^ Z I ^r Z ^r Z ^ X I ^ Z ^ Z ^r X 

+ a%a% + ololo%) . (11) 

This Hamiltonian is twisted in the sense of i.e., the 
site indices of the corresponding cluster state are ob- 
tained by the permutation (2, 3) (4, 5) ... (N - 2, N - 1) 
(cyclic notation) , for a chain of N qubits where N is even, 
see Fig. QJb,c). 

The next step in the construction of the 2D stabi- 
lizer Hamiltonian is to construct a ladder Hamiltonian 
by bonding nearest-neighbor sites on adjacent chains a 
and 6, in which all the bondings between qubits ia and 
ib are carried out simultaneously: 



ladder 



= n(-<^«0fa0&3*3a " CT 3>iVla 



CT 2a <7 4a°56 cr 5a 



alaOta<b<a) + i a b ) 



4a u 6a u 6b 

(12) 



A 2D stabilizer Hamiltonian is produced by connecting 
the above two ladder Hamiltonians with the interaction 



3 



between the two ladders. For example, when we prepare 
two ladders of length 4 such as in Fig. [ljc) and connect 
them vertically, we obtain a 4 x 4 stabilizer Hamiltonian. 

Heisenberg model- For the Heisenberg interaction, we 
can construct only a two-qubit stabilizer Hamiltonian 
(note that the same is true for the XX Z interaction). 
The basic relation is 



-i6[XYZ] l2 y je[XYZ] 12 _ 



1 



cos 2 (26>K 



+ sm z (26)a y 2 + - sin(40)(<7?<7» - a\a x 2 ) . (13) 

For the Ising and XY models, we can eliminate the 
single Pauli matrix terms leaving the interaction terms 
[see Eqs. © and ©]. However, in Eq. ([15]). if we 
set sin(2#) = or cos(2#) = 0, we also eliminate the 
o~fo~2 ~ a i a 2 term. This is because the Heisenberg inter- 
action contains terms in all three spatial directions [4] . In 
the case of two qubits, we obtain H = Q(af of — o-fa 2 ) 
from the initial Hamiltonian H s — Q(o~f — o~ 2 ) by us- 
ing Eq. (fTBf for 9 = tt/8. By applying a n- rotation, 
we obtain the two-qubit stabilizer Hamiltonian H = 



>l u 2 



Manipulation of always- on Hamiltonian - The scheme 
discussed up to now relies on switching on/off the single- 
qubit Hamiltonian Hq [see Eq. ©] and the Ising or XY 
interaction part H lni separately. There is a number of 
schemes for switching on/off interactions between qubits 
(see, e.g., [l2-14|). However, they make the system more 
complicated and require additional overhead. 

Here, we solve this problem by demonstrating how 
to extract Hq and Hi nt by using appropriate pulse se- 
quences. We illustrate the idea using the standard NMR 
Hamiltonian H nmv = £\ e^of + J2i<j J a i a j which has 
the property that [H ,H int } = 0. In this case, H and 
Hi nt can be switched on/off by using a simple pulse se- 
quence. The interaction part Hi slng can be extracted by 
using two sandwiched 7r-pulses such as exp(iT-f/i s j ng ) = 

e -i(^/2)S J ^ e i(T/2)H 1UMe i(-ir/2)Ej^ e i(-r/2)H nm . < Qn the 

other hand, two steps are required to obtain Hq. 
Let us consider a ID qubit chain. By apply- 
ing a 7r-pulse about the a;- axis (denoted by (tt) x ) 
to all the qubits on the even sites, we ob- 



tain e 



-i(*72)£< 



i(r/2)H nm r p i(7r/2)E< 



J(T/2)H aB 



giTEi^ai-i. Similarly, we obtain e lT ^ in ' 7 ^- 1 by ap- 
plying a (7r) 2; -pulse to all the qubits on the odd sites. 
Combining these two processes yields Hq- This method 
is easily generalized to 2D or 3D qubit arrays. 

If [H ,H int ] ^ 0, this NMR method cannot be used. 
Even in this case, Hq and H lnt can be extracted sepa- 
rately. The idea follows from average Hamiltonian theory 
which is based on the Baker-Campbell-Hausdorff (BCH) 
formula for the expansion of e A e B A stroboscopic 
application of the Hamiltonian designed by a series of 
short pulses can reduce or eliminate unwanted terms, if 
fir, Jt <C 1. First we extract Hi nt by setting A = ho + hi 
and B = —ho + hi in the BCH formula, where ho = Hqt 



and hi = H- lnt T. B is realized by applying a (7r)j,-pulse 
on every qubit. From the BCH formula, we obtain 



„AB 



exp(2/ii + [ho, hi 



i[/i , [h a , hi}}) 



(14) 



The exponent corresponds to a third-order expansion in 
fir for Or < 1. If we repeat this operation n times like 



e A e B e A e B , 



■e A e B 



(e 



A n B\ 



such that nfir = 7r/4, 



the fc-th term is of order [7r/(4n)] fe . Therefore, Ho is 
cancelled, and we obtain only _ff; nt in this order. When 
we apply 



e A e B e B e A 



exp(4/ii - - [hi, [ho, hi]] 



i[/io, [h ,hi]]) 



(15) 

we can eliminate the second term in Eq. (fl4]) . In the 
limit n — > oo under the condition of nflr = 7r/4, Hq is 
exactly eliminated. The extraction of Hq can be achieved 
analogously. Moreover, as shown in (lB| . if the fcth-order 
term is the first nonvanishing correction, the decay rate 
T2 of the qubit system is enlarged according to T' 2 cx 
T 2 (k + 1)1 (T 2 /t c ) k as long as T 2 > t c where t c is the time 
required for each single step [A and B in Eq. (|14[) ]. 

For the XY model, we have to switch off subsets of H lnt 
corresponding to Si and S2, as discussed after Eq. (|10[) . 
This is equivalent to choosing A and B appropriately: 
e.g., for the ID chain, A = H mii = hi e + hi Q and B = 
hi e - hi where h u = Jt([XY} 23 + [XY] i5 + [XY} 67 + ..) 
and h Xo = Jt{[XY]i 2 + [XY] 34 + [XY] 56 + ..). That 
means, B is generated by applying n pulses to qubits 
2,3,6,7,10,11,.... 

Let us consider the ID XY model with e% = in 
Eq. ((2J. The following operation can be used to obtain 
Hq- (i) Applying a (7r) a; -pulse to all the qubits on the 
even sites changes the sign of J2i<j °f °f ■ ( n ) By fur- 
ther applying a (7r) y -pulse to the same subset of qubits, 
we obtain B/t = Ei^Ki-i ~ °2i) ~ Ei<j A XY ]ij ■ 
As a result, we obtain ^X)i°2i-i- Repeating the same 
operations with the qubits on the odd sites, we obtain 
QJ2i a 2i- For the Ising model, the process (i) is not re- 
quired. For the Heisenberg model, the same procedure as 
in the case of the XY model does not eliminate the term 
JJ2i<j a i a j ■ Thus, additional similar steps are required 
to eliminate this term. 

Robustness.- Since a practical realization of these 
pulse sequences will not be free of imperfections, 
we now analyze the effect of pulse duration errors 
5. A central quantity will be the cluster-state fi- 
delity F st {r) = |(*oo...o|E/ r ( ( 5)|*oo...o)| 2 where U T (S) = 
e -jrff 3tab ( 1 r/4+5/4)^ In the Iging cas6j H stah (6) is given by 

Eq. ((5J, and an analogous equation in the XY case. 

Let us consider the ID Ising case. For O — tt/A + S/A, 
the first-order correction to Eq. §5§ reads 



H'(6) = -0 



Ki<N 



'1 + 



l) 



(16) 
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FIG. 2: Dependence of the cluster-state fidelity for r = 
n/(4Q) on random errors in the pulse duration for one- 
dimensional XY qubit arrays of length TV = 6, 8, and 10. 
The width of the Gaussian distribution is denoted by a; the 
error bars indicate the standard deviation. 



The effect of these terms is calculated from perturba- 
tion theory using the expressions for erf |^oo..o)j where 
l^oo. ..o) is the initial cluster state and a = x,y,z fltij . 
The lowest-order expression of the cluster-state fidelity 
reads F 8t (r) « [l + n 2 r 2 S 2 (N - l)/2] ~ 2 , and the correc- 
tion scales with 8 2 which is a signature of the robustness 
of our method. 

The simplest and most powerful method to further re- 
duce the effect of pulse imperfections is the symmetriza- 
tion of the pulse sequence frequently used in NMR 9]. 
We first note that Eq. is equivalent to 



P(0) 



Tlff lp , 



-Tiffin 



-Tiffin 



rffp nffipt 



P(t), (17) 



where, as before, t\ = 7r/(4J). The second half of 
this pulse sequence results in a perturbation term that 
has the opposite sign as compared to Eq. (|16l) . Ap- 
plying Eq. (|14p leads to a cancellation of the first- 
order perturbation term. If the original interval length 
r is divided into an even number n of subinter- 
vals, r = nt c , the perturbation term is replaced by 
[H s ta.b, H']iT5/(4n), and the fidelity F™ m is given by 

F£ m ps [l + n 4 (2T 2 /n) 2 5 2 (N -l)/2]~ 2 . Hence, the fi- 
delity is improved, F s ™ m > F st , if n > 2f2r. 

These pcrturbative results are complemented by ex- 
act numerical calculations of F st (r) and the gate fidelity 
F s (t) = 2~ N |Tr U}(0)U T ({6i})\ for r = tt/(40) in sys- 
tems with up to 10 qubits. Here, Si, i = 1, . . . , N — 1 
corresponds to the qubit pair + 1). In both the Ising 
and XY case, we averaged these fidelities over 2000 ran- 
dom realizations of the Si taken from a Gaussian distribu- 
tion with varying width a. The results indicate that the 
method is rather robust even outside the regime where a 
is much smaller than 7r/4 ps 0.78. For instance, the XY- 
model calculation shows that F st is bigger than 99% for 
a < 0.04, while F g can be bigger than 99% even for a as 
large as 0.07. The comparison between the two models 
shows that the robustness in the XY case (see Fig. [2J is 
somewhat better than in the Ising case. 

Discussion and Conclusion.- The stabilizer formal- 



ism relies on the multiplication of Pauli matrices and 
is widely used in quantum information theory, e.g. in 
quantum error-correcting codes (QECC) [l7|- The dis- 
tinguishing feature of the stabilizer Hamiltonians dis- 
cussed here is that their ground states are directly re- 
lated to the universal resource of measurement-based 
quantum computation. In QECC, the stabilizer formal- 
ism is used to express codewords, and the method illus- 
trated in Eq. ((6]) can be used to obtain eigenvalues for 
the syndrome measurements in the process of detecting 
errors. Moreover, relations © and ([5]) can be used to 
effectively generate codewords in solid-state qubits. For 
example, the three-qubit GHZ state (|000) ± |lll))/\/2, 
which is used for the nine-qubit code, is effective gen- 
erated by e^/^^looO), where e ^ = 

e i(7r/2)(cr^+(T|) e ±i(7r/4)crj : (Tf cr| e -i(7r/2)((T^+(T^) j g a pp nec J ^0 

the three-qubit array. The five-qubit code [17| can also 
be generated by using this method. 

To conclude, we have proposed a method to pro- 
duce 2D stabilizer Hamiltonians by using pulse sequences 
starting from natural qubit-qubit interactions. The 
ground states of these Hamiltonians are 2D cluster states 
that are universal resources of measurement-based quan- 
tum computation. We have shown that by specially tai- 
lored pulse sequences, initially prepared cluster states can 
be preserved with high fidelity. We have also shown how 
this can be implemented in the case of always-on inter- 
actions. Our work will facilitate implementations of one- 
way quantum computing. 
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project SOLID, the Swiss SNF, the NCCR Nanoscience, 
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[2 

[3] 
[4] 

[5, 

[e; 

[7 
[8 
[9 

[10 



For a review, see, e.g., H.J. Briegel, D.E. Browne, W. 
Diir, R. Raussendorf, and M. Van den Nest, Nature Phys. 
5, 19 (2009). 

H.J. Briegel and R. Raussendorf, Phys. Rev. Lett. 86, 910 
(2001); R. Raussendorf, D.E. Browne, and H.J. Briegel, 
Phys. Rev. A 68, 022312 (2003). 

M. Borhani and D. Loss, Phys. Rev. A 71, 034308 (2005). 
T. Tanamoto, Y.X. Liu, X. Hu, and F. Nori, Phys. Rev. 
Lett. 102, 100501 (2009). 

M. Van den Nest, K. Luttmer, W. Diir, and H.J. Briegel, 
Phys. Rev. A 77, 012301 (2008). 

T.C. Wei, I. Affleck, and R. Raussendorf, Phys. Rev. 
Lett. 106, 070501 (2011); 

F. Verstraete and J.I. Cirac, Phys. Rev. A 70, 060302 
(2004). 

S.D. Bartlett and T. Rudolph, Phys. Rev. A 74, 040302 
(2006). 

R.R. Ernst, G. Bodenhausen, and A. Wokaun, Principles 
of Nuclear Magnetic Resonance in One and Two Dimen- 
sions (Oxford University Press, Oxford, 1987). 
J. Dinerman and L.F. Santos, New J. Phys. 12, 055025 
(2010). 



5 



[11] R. Heule, C. Bruder, D. Burgarth, and V.M. Stojanovic 
Phys. Rev. A 82, 052333 (2010); Eur. Phys. J. D 63, 41 
(2011); V.M. Stojanovic, A. Fedorov, A. Wallraff, and C. 
Bruder, Phys. Rev. B 85, 054504 (2012). 

[12] D.V. Averin and C. Bruder, Phys. Rev. Lett. 91, 057003 
(2003). 

[13] A.O. Niskanen, Y. Nakamura, and J.S. Tsai, Phys. Rev. 

B 73, 094506 (2006). 
[14] Y.X. Liu, L.F. Wei, J.S. Tsai, and F. Nori, Phys. Rev. 



Lett. 96, 067003 (2006); M. Grajcar, Y.X. Liu, F. Nori, 
and A.M. Zagoskin, Phys. Rev. B 74, 172505 (2006). 
[15] U. Haeberlen and J.S. Waugh, Phys. Rev. 175, 453 
(1968). 

[16] H. Aschauer, W. Diir, and H.-J. Briegel, Phys. Rev. A 

71, 012319 (2005). 
[17] D. Gottesman, Ph.D. thesis, Caltech (unpublished); 

|quant-ph/9705052| 



